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Abstract. I present constraints on cosmological param- 
eters in the Ao-i^o plane from a joint analysis of gravi- 
tational lensing statistics (Helbig et al. 1999b) and the 
magnitude-redshift relation for Type la supernovae (Perl- 
mutter et al. 1999; Riess et al. 1998). I discuss reasons 
why this particular combination of tests is important and 
how the constraints can be improved in the future. The 
lensing statistics and supernova results are not inconsis- 
tent, thus it is meaningful to determine joint constraints 
on Aq and Qq by combining the results from both tests. 
The quantity measured by the lens statistics and the m-z 
relation for type la supernovae discussed here is approx- 
imately Aq — r2o. At 95% confidence, the upper limit on 
Aq — f^o from lensing statistics alone is 0.45 and from su- 
pernovae alone is in the range 0.65-0.81 (depending on 
the data set). For joint constraints, the upper limit on 
Ao — ilo is in the range 0.55-0.60 (again depending on the 
data set). For a flat universe with Aq -I- l^o = Ij this corre- 
sponds to upper limits on Aq, taking the top of the range 
from different data sets, of 0.72, 0.90 and 0.80 for lensing 
statistics alone, supernovae alone and the joint analysis, 
respectively. This is perfectly consistent with the current 
'standard cosmological model' with Aq ~ 0.7 and flo « 0.3 
and is consistent with a flat universe but, neglecting other 
cosmological tests, does not require it. 
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1. Introduction 

Recently, several papers (e.g. Ostriker & Steinhardt 1995; 
Turner 1996; Bagla et al. 1996; Krauss 1998; White 1998; 
Tegmark et al. 1998a,b; Eisenstein et al. 1998, 1999; Web- 
ster et al. 1998; Bridle et al. 1999; Efstathiou et al. 1999) 
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have pointed out the advantages of joint analyses of cos- 
mological parameters, i.e. combining the information from 
more than one cosmological test. Ideally, such tests would 
be complementary, i.e. the degeneracy in the Ao-f2o plane 
would be in orthogonal directions. However, even if this 
is not the case, indeed, even if the degeneracy is exactly 
the same, the combination of tests can tighten the con- 
straints as well as serve as a consistency check. Here, 
I discuss constraints on cosmological parameters in the 
Ao-51o plane from a joint analysis of gravitational lensing 
statistics (Helbig et al. 1999b, hereafter Paper II) and the 
magnitude-redshift relation for Type la supernovae, using 
the results of the Supernova Cosmology Project and the 
High-Z Supernova Search Team (Perlmutter et al. 1999; 
Riess ct al. 1998, hereafter SCP and HZSST, respectively). 
Although both tests are preliminary in the sense that they 
will improve with more and better observational data, the 
time is already ripe for a joint analysis, to demonstrate 
both what already can be done and how each test can be 
improved to lead to tighter joint constraints. 

The plan of this paper is as follows. In Sect. 2 I briefly 
review the basis of each of these two cosmological tests. In 
Sect. 3 I present and discuss the joint constraints. Sect. 4 
provides a summary and conclusions. 



2. Theory review 

I here use the notation of Kayser et al. (1997) with regard 
to cosmology and refer the reader there for the relevant 
definitions. In particular, f2o refers only to 'ordinary mat- 
ter' and Ao is the normalised cosmological constant, such 
that Aq + ilo = 1 for a flat universe. 

Both gravitational lensing statistics and the 
magnitude-redshift relation are 'classical' cosmologi- 
cal tests, i.e. the theoretical dependence of an observable 
quantity on redshift is compared with observations. This 
is done straightforwardly in the case of the magnitude- 
redshift relation, and in a somewhat more roundabout 
way in the maximum-likelihood analysis of gravitational 
lens statistics used here. The redshift range probed by 
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the magnitudc-rcdshift relation extends at present out 
to z 1. In the case of lensing statistics, the source 
population extends to quite large redshifts {z w 4) 
although the redshift range of significant optical depth 
is smaller. Thus, the two tests are both 'global' rather 
than 'local' cosmological tests and probe similar, though 
not identical, redshift ranges. Otherwise, the tests are 
completely independent. 

The m-z relation is concerned essentially only with 
the luminosity distance whereas lensing statistics deal 
with several different angular size distances (between ob- 
server and lens (0^), observer and source (Dg) and lens 
and source (-Dds)) (sec, e.g., Kayser et al. 1997, for a dis- 
cussion of the various cosmological distances) and the vol- 
ume; they also depend on several other 'astrophysical' pa- 
rameters (e.g. Kochanek 1996; Quast & Helbig 1999, here- 
after Paper I). 

2.1. The m-z relation for type la supernovae 

The basic idea of the m-z relation is simple: one has an 
object of known absolute magnitude M and compares it 
to the observed magnitude m. The difference or distance 
modulus is 

m - M = 5 logio D^ + K + 42.384 - 5 log^o K (1) 

where is in units of the Hubble length, K is the K- 
correction and h is the Hubble constant in units of 100 
km/s/Mpc (see, e.g., Kayser et al. 1997, for a derivation)^. 
This depends on the cosmological model since depends 
on the cosmological parameters Aq and Hq. Note that, as 
is often the case in practice, if M is known modulo h, then 
Eq. (1) does not depend on the Hubble constant at all. On 
the other hand, if M is known absolutely, this is equiva- 
lent to knowing h, assuming one has at least one object at 
low redshift (where the dependence on Aq and fio is neg- 
ligible). In any case, our knowledge (or lack of it) about 
the value of the Hubble constant Hq does not appreciably 
affect the ability of this cosmological test to measure the 
cosmological constant Aq and the density parameter ilo • 

Thus, one has a number of objects with observed mag- 
nitudes mi and a way of calculating the absolute magni- 
tudes Mi (see, e.g., SCP and HZSST for a description of 
how this is done in practice) and fits for the parameters Aq 
and Oq- If all objects are in a narrow redshift range, then 
confidence contours in the Aq-Oq plane will only allow one 
to measure approximately Aq — r^o whereas having objects 
at different redshifts breaks this degeneracy (e.g. Goobar 
& Perlmutter 1995). 



^ Note that the second occurrence of the term 'Hubble 

length' in Kayser ct al. (1997) should actually be 'Hubble 
length for /i = 1', although this is obvious from the context. 



2.2. Gravitational lensing statistics 

See Paper I and references therein for a discussion of how 
constraints on Ao-^io are derived from gravitational lensing 
statistics. Gravitational lensing statistics, at least in the 
'interesting' part of parameter space, constrain approxi- 
mately Ao — r^o (e.g. Cooray 1999). Thus the degeneracy 
is approximately the same as that of the m-z test. Thus, 
rather than reducing the allowed area of parameter space 
through orthogonal degeneracies, these two cosmological 
tests provide a consistency check on each other. Also, the 
m-z relation provides a good lower limit on Aq while lens- 
ings statistics provides an upper limit; obviously, the for- 
mer should be smaller than the latter. If this is the case, 
then the two cosmological tests arc consistent with each 
other, and it is meaningful to construct joint constraints, 
which allow a region of parameter space smaller than that 
allowed by either test alone. 

3. Data and results 

3.1. Individual results 

For the m-z test I used the results presented in SCP and 
HZSST, which have kindly been made available by the 
respective collaborations, as well as our own results from 
the analysis of JVAS, the Jodrell Bank-VLA Astromet- 
ric Survey (Paper II and references therein). Fig. 1 shows 
the likelihood ratio as a grey scale and the 68%, 90%, 95% 
and 99% confidence contours for the results from SCP and 
HZSST; Fig. 2 does the same for the results from Paper II. 
Sec these references for discussions of these results indi- 
vidually. I use the JVAS results of Paper II, rather than 
those of Paper I, since the former seem more reliable, de- 
spite the remaining uncertainties (see Paper II for a discus- 
sion). Also, using only one set of lensing statistics results, 
rather than a combination, is conservative, since the joint 
constraints are tighter than individual constraints.^ 

The basic format here is that of a probability density 
function, i.e. a relative probability as a function of Aq and 
VLq. Ideally, this would cover all values of Aq and fJo, or 
at least all for which there is a non-negligible probabil- 
ity. Alternatively, one can impose a prior constraint on 
Ao, f^o or both, such that there is a non- negligible prob- 
ability only in a comfortably small region of parameter 
space. The simplest way to do this is to use a top-hat 
function, such that the a posteriori likelihood is given by 
the a priori likelihood within some range and is exactly 
zero outside of this range. This is a conservative approach 
if the allowed range is large enough to include the cor- 
rect cosmological model in any case and also since, within 

^ Of course, if one is concerned with the consistency of the 
results, rather than in reducing the parameter space through 
joint constraints, then one should use as many results as possi- 
ble. However, it only makes sense in this context to use reliable 
results, so this is a reason to neglect the results based on optical 
gravitational lens surveys discussed in Paper I. 
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Fig. 1. The likelihood function p{D\Xo, Slo) (cf. Paper I) from Perlmutter et al. (1999) (SCP, equivalent to their Fig. 7; hereafter 
data set A) (left column), Riess et al. (1998) (HZSST) (Ami5(_B) method, equivalent to the dotted contours of their Fig. 7; 
hereafter data set B) (middle column) and Riess et al. (1998) (MLCS method, equivalent to the dotted contours of their Fig. 6; 
hereafter data set C) (right column) in the original parameter space and resolution (top row), in the parameter space used for 
the calculations in this paper but in the original resolution (middle row) and in the parameter space used for calculations in 
this paper in the resolution used for calculations in this paper (bottom row). The pixel grey level is directly proportional to the 
likelihood ratio, darker pixels reflect higher ratios. The contours mark the boundaries of the minimum 0.68, 0.90, 0.95 and 0.99 
confidence regions for the parameters Ao and fio 



the allowed range, the hkclihood depends only on the cos- 
mological tests considered and not on the priors (which 
makes for easier interpretation). 

The first three rows of Table 1 show the range of pa- 
rameter space covered by the references which are used to 
provide input data for this work. I take as the only prior 
that the likelihood is zero outside the overlap of the vari- 



ous ranges of the various cosmological tests used, as in the 
last row in Table 1. The lower limit on Jlp = is physical 
and the upper limits Oq = 2 and Ao = 2.9 are certainly 
large enough (see the discussion in Paper I on these values 
and on the use of prior information in general) . The lower 
limit on Aq comes mainly from the fact that Aq < —1 is 
strongly excluded by the m-z relation itself, although the 
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Fig. 2. The Ukehhood function p{D\\o,flo) from Paper II (hereafter data set D) in the original parameter space (left) and in 
the parameter space used for the calculations in this paper (right). (The resolution in Paper II is (due to the fact that the lens 
statistics calculations are numerically much more demanding) the worst of all the data sets considered here and is thus the one 
used for the calculations in this paper.) See Fig. 1 for a description of the plotting scheme 



Table 1. The range of Aq and fio explored by the references used here 



Reference 


range 


Ao 


resolution 


range 


no 


resolution 


Perlmutter et al. (1999) (SCP) 


-1.00 


2.98 


0.02 


0.00 


2.99 


0.01 


Riess et al. (1998) (HZSST) 


-1.00 


3.00 


0.01 


0.00 


4.00 


0.01 


Helbig et al. (1999b) (Paper II) 


-5.00 


3.00 


0.10 


0.00 


2.00 


0.10 


this work 


-1.00 


2.90 


0.10 


0.00 


2.00 


0.10 



analysis of current cosmic microwave background obser- 
vations (e.g. Lincwcavcr 1998; Macias Perez et al. 1999; 
Helbig et al. 1999a) also suggests this. In any case, an 
analysis of joint constraints from other cosmological tests, 

even excluding the m-z relation, suggests that Ao > is a 
robust result (e.g. Roos & Harun-or-Rashid 1999). 

With an individual test, likelihood contours are found 
(in all of the cases discussed here; see the discussion in 
Paper I for possible caveats when comparing the results 
of various cosmological tests as presented in the literature) 
by finding the highest contour at constant likelihood such 
that the corresponding fraction of the total likelihood is 
enclosed. Note that these contours depend not only on the 
likelihood ratio but also on the range of parameter space 
plotted. In this work, the parameter space plotted should 
be considered to have an a priori likelihood of 1, while the 
parameter space outside the plot should be considered to 
have an a priori likelihood of 0. Note that these arc 'real' 
likelihood contours, not approximations based on Ax^, 
the assumption that the probability distribution is (a 2- 
dimensional) Gaussian etc such as one often finds in the 
literature. 



It should also be noted that I take the likelihood as 

presented in the references in question. In the case of the 
lens statistics (data set V), all parameters except Aq and 
flo were held constant. In the case of the m-z relation 
(data sets A C), the results are obtained by marginalis- 
ing over the nuisance parameters (see, e.g., Paper I for a 
discussion). However, this is of no concern at the level of 
acciiracy I am concerned with here, especially since there 
are no nuisance parameters common to the m-z test and 
the lensing statistics test.^ 



^ The publicly available data from SCP are actually not the 

likelihood itself, but rather the value for each point in the pa- 
rameter space is the (normalised) sum of all likelihood values 
for all points in the parameter space which are not less than 
the value for the point in question. This format, which allows 
one to immediately plot a given confidence contour by plotting 
a contour at that level, I have converted back to the original 
probability density function. 
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Fig. 3. The 68% (top left), 90% (top right), 95% (bottom left) and 99% (bottom right) confidence contours for each of the 
data sets. The thick curves arc for the lensing statistics results (data set I>). In all plots, data set A has the contour with the 
lowest value of fio at its maximum height. Starting from this point and moving left, towards smaller values of Ao, in all plots 
one crosses first the contour of data set B then that of data set C 



3.2. Joint constraints 

The simplest thing to do when building joint constraints 
would be to multiply the corresponding probability den- 
sity functions (PDFs).^ One can then plot confidence con- 



* Of course, this must be done at the same resolution. Rather 
than interpolate the low resolution lens statistics results, I have 
reduced the resolution of the m-z results to that of the lensing 
statistics results by using only those points in the Ao-f2o plane 
which were examined in the lens statistics calculations, all of 
which were examined by both m-z tests. 



tours in the manner described above. However, it is obvi- 
ous that this is not meaningful if the PDFs are not consis- 
tent with each other, i.e. if the region of confidence for a 
'sensible' confidence level from the joint constraints does 
not overlap with the corresponding confidence level for all 
component tests. A necessary, though not sufficient, condi- 
tion for this inconsistency to exist is that the correspond- 
ing confidence contours for the individual component tests 
do not overlap. Fig. 3 shows the 68%, 90%, 95% and 99% 
confidence contours for the four data sets considered here. 
As the 90% confidence contours from all supernovae data 
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sets overlap with that of the lensing statistics, and even 
the 68% confidence contours from two of three supernovae 
data sets overlap with that of the lensing statistics, the re- 
sults from the two cosinological tests are consistent and 
one is justified in calculating joint constraints by multi- 
plying the probability distributions of the individual tests. 
Interestingly, they are most consistent at small, but not 
too small, values of Oq- The results of this are shown in 
Fig. 4. 

Note that if there is some offset between the allowed 
regions from each individual test, as is the case here, then 
certain aspects of the joint constraints, such as in this case 
the upper limit on Aq, will not necessarily be tighter than 
the corresponding aspect from each individual test. The 
joint constraints are nevertheless better in that the allowed 
region is smaller and that this allowed region should con- 
tain the correct cosmological model, assuming of course 
that the results of the individual tests are correct as far 
as they go. The tests are very different in nature and one 
should not expect the form of the probability density func- 
tion to be the same in each case. In particular, as lensing 
statistics is especially sensitive to a large cosmological con- 
stant, the gradient in this area of parameter space is quite 
steep, thus it is not surprising that the lensing statistics 
upper limit on Aq is tighter. The fact that the confidence 
contours from the individual tests overlap shows that the 
tests are not inconsistent, and of course the allowed region 
from the joint constraints, which is consistent with each 
individual test, is approximated by this overlap. 

Since the two m-z results are not completely indepen- 
dent, the question of the consistency of or joint constraints 
from the two supernovae data sets will not be discussed in 
this paper. Rather, the question is the consistency of and 
joint constraints from each of these data sets individually 
with the lensing statistics constraints. 

Fig. 4 is the main conclusion of this paper. Although 
lensing statistics and the m-z relation individually allow, 
for appropriate values of Aq, rather large values of Qq, 
the joint constraints clearly indicate a lower fig, in ac- 



cordance with observational evidence which measures flo 
more 'directly' (see the discussion in Paper I). Compared 
to the supernovae results, the allowed region of parame- 
ter space is shifted somewhat towards lower values of f2o 
in the joint constraints. Although the actual best-fit value 
should not be taken too seriously, it is comfortably close to 
the current 'standard cosmological model' with Aq » 0.7 
and fio ~ 0.3 

The quantity measured by both the lens statistics and 
the m-z relation for type la supernovae discussed here is 
approximately Aq — JIq. Table 2 shows the 95% confidence 
ranges for Aq — Op allowed by each of the four data sets in- 
dividually and by the joint constraints of data set V with 
data sets A, B and C. At 95% confidence, the upper limit 
on Aq — rio from lensing statistics alone is 0.45 and from 
supernovae alone is in the range 0.65-0.81 (depending on 
the data set).^ For joint constraints, the upper limit on 
Aq — ^0 is in the range 0.55-0.60 (again depending on the 
data set). For a flat universe with Aq -|- fio = 1, this corre- 
sponds to upper limits on Aq, taking the top of the range 
from different data sets, of 0.72, 0.90 and 0.80 for lens- 
ing statistics alone, siipernovac alone and the joint anal- 
ysis, respectively. Again, this is perfectly consistent with 
the current 'standard cosmological model' with Aq w 0.7 
and rig ~ 0.3 (e.g. Roos & Harun-or-Rashid 1999; Turner 
1999) and is consistent with a fiat universe but, neglecting 
other cosmological tests, does not require it. 

3.3. Systematic errors 

As far as the m-z relation for type la supernovae goes, 
various possible sources of systematic errors have been 
discussed in detail by SCP and HZSST. See particularly 
Fig. 5 in SCP. Basically, there is no evidence that the 
purported effects could significantly bias the results or, as 

^ Note that this is not the same as that quoted in Paper II; 
this is because, as discussed above, the range of parameter 
space examined or, equivalently (in this case), the prior is dif- 
ferent. 
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Table 2. 95% confidence ranges for Ao — f^o allowed by each of the four data sets individually as well as various joint constraints 



Reference 


lower limit 


upper limit 


Perlmutter et al. (1999) (SCP) (data set A) 


-0.05 


+0.65 


Riess et al. (1998) (HZSST, Ami6(S)) (data set B) 


+0.30 


+0.81 


Riess et al. (1998) (HZSST, MLCS) (data set C) 


-0.12 


+0.79 


Hclbig et al. (1999b) (Paper II) (data set V) 


-1.90 


+0.45 


A + V 


-0.20 


+0.55 


B + V 


±0.00 


+0.55 


C + V 


-0.25 


0.60 



in the case of grey dust, while a modest amount cannot 
be ruled out, it seems physically rather implausible (but 
see, however, Aguirre 1999a, b; Aguirre fc Haiman 1999). 
(Note that Falco ct al. (1999) find no evidence for grey 
dust at optical wavelengths based on studies of extinction 
in gravitational lens galaxies.) It is interesting to note that 
while one can invoke grey dust to explain the dimming of 
supernovae relative to the expectation in for example a 
dust-free Aq = model, instead of invoking for example a 
low-density model with a positive cosmological constant, 
this degeneracy can be broken by observing supernovae at 
higher redshift (e.g. Aguirre 1999a) than has been done 
up until now: the m-z relation as a function of Ao and 
Oo is exactly known, so the larger the range in redshift 
for which the m-z relation is observed, the more ad hoc 
alternative explanations become, provided of course that 
there is a cosmological model (which is not ruled out on 
other grounds) which provides an acceptable fit to the 
data. 

Grey dust is also something which can effect gravita- 
tional lensing statistics based on optical samples, although 
the effects are not so straightforward. On the one hand, 
if the grey dust is concentrated in (lensing) galaxies, this 
could lead to lens systems being missed in the siirvcy. To 
first order, this would lead to an imderestimate of the 
optical depth and thus of the value of Aq. On the other 
hand, again if dust is concentrated in (lensing) galaxies, 
the sources in the identified lens systems can suffer from 
extinction, which, depending on the details of the luminos- 
ity function, could lead to a wrong estimate of the mag- 
nification bias. In the 'normal' case of a fiattening of the 
luminosity function for fainter objects, this will lead to 
an underestimate of the amplification bias and hence an 
overestimate of the optical depth and thus the value of 
Ao (Falco et al. 1999). Radio surveys of course are not af- 
fected by dust, so in principle one could detect grey dust 
through a systematic difference in the results from optical 
and radio surveys. In practice, however, the presence of 
other systematic effects makes such a detailed comparison 
impractical. 

Radio surveys for gravitational lenses offer many ad- 
vantages over optical surveys (see the discussion in Pa- 
per II). However, at present, the main source of uncer- 
tainty, lack of knowledge about the source population, 
makes them worse than optical surveys in this respect. 



For the calculation of the amplification bias, one needs 
to know, at a given redshift, the luminosity function.^ 
On the other hand, much information is gained from the 
sources in a survey which are not lensed (see the discus- 
sion in Paper I); to interpret this, one needs to know, at a 
given flux density, the redshift distribution of the sources. 
Of course, these two things — the redshift-dependent lu- 
minosity function and the flux-density dependent redshift 
distribution — are different sides of the same coin. 

The number counts of the Cosmic Lens All-Sky Sur- 
vey, of which the Jodrell Bank-VLA Astrometric Survey, 
the results of which are used here, is a subset, suggest that 
amplification bias is not a big effect; hence the systematic 
error from lack of knowledge of the luminosity function is 
probably small, although it is conceivable that the number 
counts (integrated over redshift, as in general the CLASS 
sources have unknown redshifts) of CLASS are not rep- 
resentative of the luminosity function at all redshifts. As 
the lensed sources are generally at higher redshifts, the 
luminosity function might be different here and thus the 
true amplification bias different from that which was used 
in the JVAS analysis of Paper II (where it was assumed 
that the CLASS number counts are representative of all 
redshifts). 

In Paper II, it was also assumed that the redshift dis- 
tribution of JVAS is equal to that of CJF, independent 
of flux-density. There is some preliminary evidence that, 
as one moves toward lower flux-density levels, the typical 
redshift of flat-spectrum radio sources decreases. If this 
is the case, then our JVAS analysis will have underesti- 
mated the value of Aq, as a higher value of Aq (all other 
things being equal) is needed to achieve the same opti- 
cal depth for a low-redshift source than is needed for a 
high-redshift source. Although the results from the m-z 
relation for type la supernovae and gravitational lensing 
statistics are not inconsistent, and although, due to the 
different dependence on the cosmological parameters, the 
lower limit on Aq will always be stronger from the former 
and the upper limit from the latter, if it does turn out 
to be true that the CJF redshift distribution is systemat- 

® Due to the amplification of the gravitational lens effect, 
lensed sources near the lower flux-density limit of the survey 
will have an unlensed flux density lower than this, so the lumi- 
nosity function thus needs to be known down to a flux-density 
level a factor of several below that of the survey. 
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ically higher than the true redshift distribution of JVAS, 
then the results from the m-z relation for type la super- 
novae and gravitational lensing statistics will become even 
more consistent. 

4. Summary and conclusions 

I have presented the first detailed analysis of joint con- 
straints between gravitational lensing statistics and the 

m-z relation for type la supernovae, making use of data 
from Helbig et al. (1999b), Perlmutter et al. (1999) and 
Riess et al. (1998), presenting the individual results and 
the new joint constraints in a uniform way. The two tests 
are not inconsistent, the joint constraints are tighter than 
those from either test individually and provide additional 
evidence in favour of the current 'standard cosmological 
model' with Aq « 0.7 and fio ~ 0.3, although (neglecting 
constraints from other sources such as the CMB) a reason- 
able range of other cosmological models is not excluded. 

In the near future, gravitational lensing statistics from 
CLASS, the Cosmic Lens All-Sky Survey (Myers et al. 
1999) should reduce both the random and systematic er- 
rors. Should the results from lensing statistics and the m- 
z relation for type la supernovae remain consistent, this 
should reduce the allowed parameter space even further. 
We are truly entering an era of precision cosmology, where 
the overlap of the allowed regions of parameter space from 
many different and independent cosmological tests is very 
small but not zero. 

The data for the figures shown in this paper are avail- 
able at 

http : //mult i vac . jb . man .ac.uk: 8000/ceres 
/ data_f rom_papers/ snl ens /snl ens .html 

or 

http : //gladia . astro . rug . nl : 8000/ceres 
/data_from_papers/snlens/snl ens .html 
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